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. DIFFERENTES APPROCHES DE L' ANALYSE SPECTRALE 


RESUME : Depuis quelques annees, plusieuvs approches du 

problime de la recherche de la densite spectrale 
de puissance d'une fonction aleatoire a. partir d'une estima- 
tion de la fonction d 'autocoi‘relation ont ete itudi^es. 

Nous presentons id une etude comparative de ces 
differentes methodes, en mettant en evidence les prinoipes 
SUV lesquels elles reposent et les hypotheses qu'elles im- 
p liquent. 


Nous donnons egalement quelques indications sur I' op 
timisation de la "longueur" de la fonction de correlation 
estimSe et nous terminons par un exemple d' application des 
diffdrentes mdthodes prisentees. 


ABSTRACT : DIFFERENT APPROACHES OF SPECTRAL ANALYSIS . 

In the last years, several approaches of the problem 
of the calculation of the spectral power density of a random 
function from an estimate of the autocorrelation function have 
been studied. 

We present in this paper a comparative study of 
these different methods and point out the principles on which 
they are based and the hypothesis implied. 

We give also some indications on the optimization 
of the "length" of the estimated correlation function and 
finally we give an example of application of the different 
methods discussed in this paper. 
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1. POSITION DU PROBLEME . 

L 'observation est une realisation particuli^re 
d'une suite aleatoire : 





X(l) . . .X(M) . . . = { X(n) } 


dont la fonction de correlation estimee est : 

. f^(k) |kl s N 

Nous envisagerons successivement les deux questions suivantes : 
Probleme 1 : 

rck) j |k| N etant supposee connue sans erreurs , trouver 
la function spectrale de {x} ; PCX) telle que 

p+TT , 

eiXk 
-TT 

Pr'obleme . 2 : 

On ne dispose que d ' estimations ^j^Ck) de C (k) s N] 

calculees S partir d'une suite de longueur M. Trouver le N 
optimal. 


I' (k) = 

^ /2iF 


5. EXEMPLES DE PROBLEMES PRATIQUES CONDUISANT A GETTE SITUATION . 
2 . 1 Analyse spectrale temporelle . 

Tx(k) = r^^CkT^) 

fonction ' de correlation echantillonnee.de la fonction aleatoire 
XCt). 
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■ La densite spectrale de puissance de la fonction 
aleatoire X(-t) ' est : , 




Yj^(v) = TF[Tj^(t)] 


Si Yy(v) = 0 pour v ^ 2T~ ^ spectre borne) 


P(A) = — ^ YyC ) 

T /2? ^ 

Ij 


pour - 7T < X < +'ir 


2 . 2 Antenne formee de multioapteurs . 

Sbit une antenne lineaire formee de N capteurs omnidirec- 
tionnels S^, Sj^. Les sources a 1' inf ini emettent un 

signal monochromatique (frequence v^). L'amplitude du signal 
emis par les sources est une fonction aleatoire de 1' angle 
0 : F(0). 



On suppose que le champ de sources est a correlation spatiale 
microscopique , soit : 

; E[F(0) F*(0')] = H(0) 6(0-0') 
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H( 0 ) : intensite dans la direction 0 . Alors 


G(3) = H Arcsin 77^ 

2 air 
J 


/T^a y'l-(M) ‘ 
2ira 


ou : 


longueur d’onde 
distance entre capteurs 


2 a sin 0 
X 


et r(k) = ECS. S. , ] 

1 i-k-^ 


_ S. 


^ : sortie de capteur i 


r(k) = interco^^^elation des signaux regus sur deux 
capteurs distants de Ka 


sont relics par 


G(3) e^^^ d3 


a - A 
2 


la relation (2) se ramene a (1) 


a < A 
a < 2 


on se trouve dans une situation analogue au sur- 
echantillonnage . Lorsque 0 varie de -it/2 a 


+tt/2 , 3 varie de — jp-rr 


2a .2a^.x 


suffit, pour se ramener a (1), de prolonger 

r\ _ r\ _ 


G(3) par 0 sur C-ir, ~) UC^tt , tt) 


a > ^ 


lorsque a varie de a , 0 varie de 

-^ir (<-ir) a (>+tt) 


On peut alors se ramener a (1) seulement dans le 
cas oii G( 3 ) = 0 pour l 3 l > 'U' , (les sources soi 
confinees dans une region de I'espace). 
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(1) Lorsque 3 var-te de -tt/ 2. d +ir/2j t'antenne balaie tout 
Z ’espaae 


(2) Lorsque 3 varie de -32 <5 t -^2 » t'antenne balaie tout 
I 'espaae 

(3) Lorsque 3 varie de -T!/2 d ■t‘n/2 , t’antenne balaie une 
partie de I'espaae. 

FIGURE 2 


Lorsque les conditions definies ci-dessus sont 
remplies, le passage de 1 ' observation CTCk)] a I'inconnue 
[H(6) ou G(3)3 se ramene au probleme d'analyse spectrale 
defini en 1. 


3-. RECHERCHE DE LA d.s.p. A PARTIR D'UNE SERIE LIMITEE DE 
VALEURS EXACTES DE LA FONCTION DE CORRELATION. 

r (k) connu sur -N ^ k ^ N 

Soit : z = e^^ 
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Q(z) = PC-i Log z] : transformee en z de 


r^(k) = 


i/2¥ 




Q(z) z.^ ^ dz 


cercle unite 




Q(z) = 


/2it 


+ N , r- /-I ^ 

I * I ^ 

-N lkl>N 


Q(z)=Q (z)+Q. (z) 

^ ^connu ^inc. 


- Q (z) est construit a partir des valeurs connues de TCk) 

connu 

- Q. (z) depend de la partie inconnue de F(k) 

me . . 

Les differentes methodes consistent a choisir, selon 

certains "criteres", une valeur pour Q,. (z). Le fait de fixer 

^ inc 

Q^j^^(z) revient a "pif'olonger" F(k) au dela de k=N. 


3.1’ Transformee de FOURIER. 


On choisit : Q. (z) = 0 

^inc. 


Alors : T ^^(k) = Tj^(k) 

= 0 


|k| s N 
lk| > N 




P C X) = - 1 - ^ r ck) e 

^ -N " 


iXk 


n'est pas s 0 \/ X dans tous les cas. 


* Solution : 


Fonction d ' apodisation (BARTLETT . . . ) , mais alors 
on ne conserve pas les valeurs connues de Fj^Ck) i 


3 . 2 Maximum d'entropie ou methode autoregressive . 

On impose : . 

^ME connu “ ^connu 

2) "L '‘entropie" H du processus prolonge est maximum. 

3.21 • 

Considerons une suite aleatoire de longueur P + 1 


{X}p = {X(0) ; . . . , X(P)} 

{X}p est une variable aleatoire a P + 1 dimensions''. Admettons 
que : 


ECX(i)] =0 Vi 

ECX(i) X(j)] = r. . 

1 J 

{X} gaussienne 


L'entropie. de {X}p est : 

~ 2 11^11 ~ 2 1 ir ^ 11 determinant de F 


( r(-F) ... r(-p) 

r(P) r(0) 

X. valeurs propres de T (s 0). 

P inf ini, on definit la "densite d'en- 


et la ’’densite d'entropie" d'une suite «> en faisant tendre 


Pour passer au cas 
tropie"i 


H- 


H 


■P. 


P P+1 



P vers I'infini 


H = lim H. 
P-»-oo ^ 


H = 


lim 

P-+-00 


i Log 

p + 1 


Un theoreme de SZEGO nous indique que : 


t Log X . +TT 

A ^ A f 


lim 

n-»-oo P 


Log[P(X)] dX 


-TT 


d'oli I'entropie de la suite {X} 

f+% 


H a 


LogCP(X)3dX 


-7T 


3.22 RoiIi§tgue_sur_le_f ondement_du_crit|re_du_maximum_d_J_entropie 

La "validite" de ce critere s'appuie sur le fait 
que le prolongement de la fonction de correlation ainsi obtenu 
correspond au processus "le plus desordonne possible" compatible 
avec les donnees et I'on a appris en thermodynamique que les 
systemes tendent naturellement vers cet etat de desordre maximum. 

On peut aussi, plus simplement, considerer la minimi- 
sation de H comme une contrainte raisonnable (parmi d'autres) 
imposee I la partie inconnue de la fonction de correlation. 
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3.23 Mise_en_oeuyre_du_critere_du_inaximum__d_J^entropie . 


;+TT 


H = 


Log Pj^g(A) dA = i 


-IT 


Log[Q„j,(z)3 ^ 

4> 


on veut maximiser H avec la condition : 




'‘ME connu 


(z) = Q 

connu 


= Qf^£ connu 


(z) + I qj 

lkl>N ^ 


on impose done 


i. 


3H 

9qk 


= 0 


k[>N 


Soit 


z 


dz-.^: Q 


+N 


soit : 


Qme^ 


I q' 


n 


-N 


n 


k >N 


La transformee en z de 


Tr~ r T'i degre N, 


qui, combine avec 


^ME connu^^^ ~ ^connu^^^ 


fixe Qj^ig^z) • 


3.24 ^ • 

Soit Gj^ le filtre lineaire A.R. predisant 

avec une erreur quadrat ique moyenne minimum lorsque Xq, 

X. ^ sont connus et soit F.^ le filtre J-G„ 

- N+1 N N 


F^{{X}]i = X. - S„[{X}]. = X. - Xj 



1 ) 


7 . 


N 


Soit : F^Cz) = I (B^, 


• I ' invepsei <16 la'-'tp.ans forme e en 


-a 


de la reponse- impulsionnelle du filtre (r.i. de Fj^ ) 


N 




1/ 


\ 




est donnee par F . ^ = 


f A 

0 

\o J 


ou 


r = 


r(0) 


rcN) 


r(-N)\ 


rc 0) 


) 


et + ^ I'erreur quadratique moyenne de prediction. 


•■f-a- 

‘ i!' 

I 

V. 


2/ 


« , N N+1 

0 (z) = 

ME /7i? F.,(e) F„*(i,) 


N 


N 'z’ 


t-' 

'1 


Ces resultats nous donnent : 


- une methode de calcul de ou 

de YULE-WALKER) ; 


(algor ithme 


- une interpretation "physique” 
maximum d'entropie est la d.s 
engendree par le filtre Fj^ 
puissance Pn+i//^ 


de E^_(l) la d.s.p. du 
.p. de la sequence aleatoire 
excite par un bruit blanc de 


1 / 

* Cette identiti resulj;e^du fait que ^Q^^(z) et • 

N ( ^ sont de degre N ' et que 


=: 




r'N ^ 1" 


^ME^^^ " ^connu^^^ 'd I'ovdve N . 


10 , 


On peut ramener un probleme d' analyse spectrale a un probleme 
d ' identification d'un filtre lineaire. 


3.3 Method© de PISARENKO. 


On a vu que : 




+N 


-N 


q * z 
^ n 


n 


q' = q' 

^ n ^ -n 


Soit 


N 


1 


In I 

n = l 


* s 
z ) 
n 


avec 


z i 1 
n ' 


V n 


Avec cette ecriture il vient pour 




i/2v 


^ . dz k 


N 


= I q^jlzjjl cos[k(J)„] 


N'“N^ 


n 


en posant 


z = z e 
n ' n' 


i(J)n 


avec |zj^I ^ 1 


Tj^gCk) est done la sonune de N cosinusoides amorties 


PISARENKO a montre que I'on peut ecrire : 

N/2 


r^(k) = Uq 6 (k) + 


7a • cos • kd) . N pair 
1 ^ 1 


(N--l)/2 
c+ I a. cos k((). 


N impair 
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n 


avec lia condition : 


^p(k) = 


iKl s N 



Interpretation : PISARENKO interprete comine forme par 

- un bruit blanc CviQ6(k)] 

- N "raies pures" de "frequence'- 

En termes d'antenne : un bruit blanc omnidirectionnel et N 
sources ponctuelles. 


• Determination de y q (puissance du bruit blanc) et des , (|) ^ . 

Prenons le cas d'une antenne et realisons la "formation de voies" 
decrite ci-dessous : 


FIGURE 3 



^ go ^ 


g = 


La puissance en sortie est 


E S 


i-r-S 


/ r(0 ) - r(-N)\ 


\r(N) - r(p) / 


\%/ 
= 1 


• Cherchons le qui .’'.elimine les N sources ponctuelles". 

Pour ^ = gjjj ■ . EC[S|^]min = y^ . , 

est le vecteur propre associe a la valeur propre minimum 
de r . (supposee s 0) . 
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1 



I 


1 





Mg = valeur propre minimale de £ .. 


N 


ikXj _ 


Enfin : I = 0 puisque les sources sont eliminees, 


Les z. = e 

: 


iXj 


sont done les N racines de : 




3.4 Methodes hybdrides (M.E. - PIS) 

PISARENKO donne le resultat suivant : 

• Soit £. associee a la sequence {X(n)} . Soit \i tel 

que 0 < y < lig (v.p. min. de C )• Soit ; = £ - yl ' 

( £^ correspond .a une sequence dont on a retire une partie 
. du bruit blanc). 


Revenant a notre interpretation de 3.3, on a le comporternent 
suivant : aux N cosinusoides amorties de correspondent 

racines de module < 1 .. Lorsque y Mg » les N racines 
tendent vers N points de cercle unite (Zp^g). 


FIGURE 4 



z complexe 


PIS 


• D'oii les methodes "hybrides” : on remplace . £ par £ (y<yo) 


et I'on calcule Pw„ (X) associee a £ , 

MEy y , 



C E P H A G 




Comportement d'une raie : 

Le'fait de retirer une partie du 
bruit blanc" entraine une "ampli- 
fication" de plus en plus poussee 
des raies spectrales 

(On tend vers des impulsions de 
DIRAC lorsque y -»■ Pq). 



4. QUELQUES INDICATIONS SUR LES ERREURS D’ ESTIMATION ET LA 
RECHERCHE D'UNE "LONGUEUR" N OPTIMALE. 


On dispose de 
(x(n) j 0 < n < M} 


■ • 


•et on estime la fonction de correlation par : 

{Tj^Ck) ; 0 < k ^ N} 

dans les differentes situations decrites ci-dessus, nous 
allons chercher a definir et a calculer une valeur optimale 
pour N , longueur de la fonction de correlation estimee. 


4.1 Transformee de FOURIER . 

Avec I'hypothese que le support de la fonction de cor- 
relation est de I'ordre de N [ | T^(k) | <<F^(0) si |k[ > N] 
et en appelant Py(A) la vraie valeur de la d.s.p., la 
variance de I'estimee P^p(X) de la d.s.p. est : 

VarCP^pU)] = ^ p2(X) ' , 

IVerreur quadratique totale est : 
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j . 



E QT(X) = (Biais)^ + Variance 






Soit 




E- QT(X) = CPj^p^(X) - Py(X)]2. + var Pj^p^CX) 


Prenons comme mesure de I'erreur totale : 


( +7T 


Et(N> = 


EQT(X) dX 


-TT 


/:+^ 


= 2 ? 


CPojt(X) - Pv(X)]“ dX - ^ 


-IT 


+U 


Py^ (X) dX 


-ir 


qui s'ecrit, en utilisant la relation de PARSEVAL : 


+00 

E^(N) = I |rj,(k)|2 + ^ Ilr^ool" 
|k|>N 

= 1 2 

1 decroit lorsque N croit (fig 6) 

2 croit lorsque N croit (fig 6) 

d'ou 1 'apparition d'un N optimal. 



7?Tnrir)i? 2 ? 
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Exemple : 


Soit : r(k) = ql 


lk]>N ^ 


I lr(k)l 


Soit 


E^(N) 




N ^ est obtenu pour : 
opt 


dE^(N) 


soit 


log k 


opt 2 log q 


1+q^ 

^ = 2M(- Log 


q = .0,8 


M = 100 
M = 2 0 0 
M = 50 


N . = 7,4 : 8 
opt 

Nopt = 8,96: 9 
^pt = ^ ^ 


/2 

q = 0,707 = ^ M = 100 

M = 50 


Nopt = 5,53: 6 
Nopt =4,53: 5 


4.2 MEM 


AKAIKE a introduit le F.P.E. (final prediction error). 


Soit un filtre A.R. d'ordre N estimant Y, a partir d'une 
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~ n - 


estimation du filtre predicteur issue de T 


N 


FPE = E[ - Y^)^l 


N 


" I ®M,m 


m 


filtre de prediction 

M,m 






: erreur d' estimation du 


'M,m 


filtre de prediction, 


3m ■ = Bx/ + ^3m 
M,m M,m M,m 


En supposant A3j^j et Yj^ independants il vient : 


N 


N 


F.PE = E[(Y^ - I + EC( I 4B„ „ V^„)M 


M ,m km 


^ V_ 


1 erreur de prediction 

2 • erreur d’ estimation du filtre A.R. 



AKAIKE 

(FPE)„ 

N 



propose : 

_ M + CN+1) ; 2 
M - (N-1) 


estimation de 1' erreur de 
prediction. 


5. CONCLUSION. 


Dans cet expose nous avons decrit,en essayant de 
montrer les liens qui les unissent, un certain nombre de 
methodes d 'analyse spectrale. Dans ce sujet, tres vaste, 
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nous nous sommes limites aux methodes basees sur le theoreme 
de WIENER-KINTCHINE pour lesquelles la densite spectrale de 
puissance est obtenue a partir d'une estimation de la fonction 
de correlation. 

Nous persons avoir mis en evidence 1.' aspect relatif 
des "qualites" des differentes methodes et, en particulier, 
nous esperons avoir conti'ibue a demystifier la "vertu supreme" 
consistant a faire "sortir des raies pures". II est clair, 
en effet, que ce simple critere ne ,peut Servir de jugement 
object if lorsque I'on sait que la methode de PISARENKO conduit 
toujours a des raies pures. 

Quant a la question de savoir quelle est la meilleure 
methode, notre reponse sera que cela depend des connaissances autre 
que celles issues de I'estimation de la fonction de correlation 
dont dispose 1 * experimentateur . Si le signal etudie est 
eff ectivement issu du filtrage d'un bruit blanc, la methode 
du maximum d'entropie est la mieux adaptee; si le signal 
est forme de raies pures, la methode de PISARENKO s’ impose. 

Ceci nous conduit, en terminant , a insister sur 1 ' importance 
de la connaissance d'un modele du systeme physique etudie. 




/ 
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APPENDICE : EXEMPLE D ’APPLICATION 


■•O 




A.l 


amorti : 


rCk) connu sur -3, +3 et forme d'un cosinus 


rco) = 1 rci) = q 5 rc2) = o 


/2 ,3 
2 ^ 


rCk) = q^ cos[(-^)it] 


Calculons les differentes dsp definies ci-dessus. 


A. 11 Transformee de FOURIER (TF) 

' ■ ■ ^ o f 

P „(A) = -[1 + q/2 cos X - q /2 cos 3X] 

/2lf 

On constate (fig. 8 ) que cette dsp n’est pas a 0 

Calculons la dsp obtenue avec une fonction de ponderation de 
BARTLETT : 

. . Fg(k) = (1 - M 


, ^BARTLETT ^ ^ ^ 


1 

/2¥ 


Cl + 


3/2 

4 


q 


cos X - 


Q ^/2 

4 


cos 3X] 


A. 12 Maximum d'entropie : 
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r . 3 = 


/P 4 V 

0 

W 


e = 


/ 1 
3. 
3' 


Posant : FCz) : 1 

il vient 


+ q(l+^-) 

_ ^ 

2 

3 /"2 
- ^ — 

2 4 6 

1 _ - a_ _ a_ 

^ 2 4 4 

.2 


3 


32 z 


S3Z- 


- 3 _ - 


MEMCX}, 


- a_ 

4 4 


2- ^ 4 6 


6 • 


■ _3 


-/2q(l+q^)cos X+(q^-. ^ - ^)cos 2X+ cos 3X3 


q,. 4 4 


4 4 


A. 13 Methode de PISARENKO. 


On reprend F defini en 3.52 


r-p I 


/i-y 

/2 

2 

0 

^ n 
■ ~ ^ 

! r- 

q/2 

2 

1-M 

^ n 
2 ^ 

0 

0 

1 

q/2 

2 

1-y 

2 

\ n3 

0 

2 

1-U 


D*ou on deduit .les mineurs de la premiere, ligne 
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^11 ^ - q^[l-p3 

^12 = ^ qC (l-M)^ ■ ^ 

.2_^4 

A^3 = (1-y) [3-^] 

Ai 4 = - ^ q^[(l-p)2 - I - V 

et 1' equation caracteristique : 

2 6 4 

2 3q + q q 2,2 

X - 2 x + ^(l+q) = 0 

X = (1-u)^ 

Soit : 

2 

X = ^ [3 + + / B] 

B = 5 - 8q^ + 2q^ + q® 
et les 4 valeurs propres. 



Nous avons mene le calcul jusqu’au bout dans un cas 
particulier (of A2 ) . 


A . 14''' Methodes hybrides . 

Soit 0 < y < 
= r - PI 


p . 
min 
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5 k/- 


/2 3 


^ — 


/2 3 


1 1 2(l-y) 2[(l-y)^ - ] 


= q/? pC l-y)^ - q^(l+q^) 


l^(l-q^) 


2[(l-y)2 - q^] 


Posant 


Enf in 


. _ q^/2 . 2(l-y)^ - (l+g^! 

3 2 (l-y) ' 2 C(l-y)^ - q^] 


X = ( 1 -y) 


_ 4 x^-( 6 q^ + 2q^)x + q‘^(l+q^)^ 


4t^[x - q"^] 


^y 


/TttCA + B cos X + C cos 2 X + D cos 3 X] 


9 9 9 

A = 1 il- + Y2 + Y3 


B = 2 tY^ + Y^ Y2 + Y2 Y33 


C = 2CY2 + Y^ Y3] 


D = 2 Y. 
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A.2 Exemples numeriques : • 

Dans le cas decrit en A1 nous avons f ait : ^ 

1 

q = — 

Z? 

Soit 

r ( 0 ) = 1 
r(i) =0,5 
r( 2 ) =0 
r(3) =-0,25 

*•'4 

alors 

^ C 1 + cos X- 0,5 cos 3X] 

./2if 

+ 0,75 cos X = 0,125 cos 3A] 

/27T 

p f X ) • -*■' 0 ,65625 

mem' ^ ^ 1,46875 - 1,5 cos X - .0., 34375 cos 2 X + 0,25 cos X 

Ppjg(X) = 0,25 + 0,25 6(0 ) + 0,5 6(X-0,33) 

_! ^ 

A + B cos X + C cos 2X + D cos 3X 

1 

/2tt 


Ptfb^^^ = 


^MEMiy^^^ 
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Po “ 0 » 2 5 



V = 0,225 
90 % de Po 

p = 0,2375 
95 % de Pf) 

p = 0,2475 
99 % de Po 

^4 

0,171 141 

0,100 94 

0,2475 

A 

4,8651 

6,58320 

9,08697 


B 

- 7,677 

- 10,51651 

- 14,54957 

C 

3,8699 

5,32894 

7,30752 j 

S ■ 

D 

- 0,9577 

- 1,35804 

/ 

- 1,84268 


Les d.s.p. obtenues sont presentees sur les figures 
suivantes. ’ 
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